Two-channel model of photoassociation in the vicinity of a Feshbach resonance 
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We derive the two-channel (TC) description of the photoassociation (PA) process in the presence 
of a magnetic Feshbach resonance and compare to full coupled multi-channel calculations for the 
scattering of 6 Li- S7 Rb. Previously derived results [P. Pellegrini et ai, Phys. Rev. Lett. 101, 
053201 (2008)] are corrected. The PA process is shown to be fully described by two parameters: 
the maximal transition rate and the point of vanishing transition rate. The TC approximation 
reproduces excellently the PA transition rates of the full multi-channel calculation and reveals, e.g., 
that the enhancement of the rate at a resonance is directly connected to the position of vanishing 
rate. For the description of two independent resonances it was found that only three parameters 
completely characterize the PA process. 



I. INTRODUCTION 

The phenomenon of a magnetic Feshbach resonance 
(MFR) is widely used for the manipulation of systems of 
ultracold atoms. One field of interest is the combination 
of the photoassociation (PA) of molecules with MFRs. It 
has been shown, both theoretically and experimentally, 
that the PA transfer rate can be significantly increased 
in the vicinity of an MFR [3, 0, 0, 0, 0] . This leads to the 
prospect of creating a large number of ultracold molecules 
out of a sample of ultracold atoms. These molecules are 
of great interest for applications in quantum information 
processing 0, 0], the exploration of lattices of dipolar 
molecules [H , or ultracold chemical reactions d, llOf . 

For all PA schemes that exploit the enhancement of 
PA at an MFR, the understanding of the interplay be- 
tween both processes is important. Here, we seek to de- 
scribe the process by a two-channel (TC) approximation 
[TTI ] . In Q this approximation has recently been used to 
predict the behavior of the PA transition rate as a func- 
tion of the scattering length. We review this approach 
and find a simplified expression with only two instead 
of three free parameters. These two parameters which 
can be, e.g., the maximal transition rate and the po- 
sition of the minimal transition rate can be obtained ei- 
ther from multi-channel (MC) calculations or from exper- 
imental observations. They can serve as a classification 
of transition processes and reveal a universal dependence 
of the enhancement of the transition rate on the position 
of vanishing transition rate. 

The Hamiltonian of relative motion for two colliding 
ground-state alkali-metal atoms is given by [l2j 

2 

H = T^ + ^(Vf + V?) + V int (1) 

where T M is the kinetic energy and fj, is the reduced mass. 
The hyperfine operator V^ f — -^tSj ■ ij and the Zeeman 
operator = (j e Sj — Jnij) ' B in the presence of a mag- 
netic field B depend on the electronic spin Sj , the nuclear 
spin ij, the hyperfine constant a J hf of atom j — 1,2, and 



on the nuclear and electronic gyromagnetic factors 7„ 
and 7 e . The central interaction 

OR) = V (R)P + Vi(R)Pi (2) 

is a combination of singlet and triplet Born-Oppcnheimer 
potentials Vq(R) and V\(R) where Po and Pi project re- 
spectively on the singlet and triplet components of the 
scattering wave function. 

For low collision energies the eigenfunctions of Hamil- 
tonian ((T|) may be written as a superposition 

|*>=5>a(i*)l«> (3) 

a 

of s-wave functions in the atomic basis |a) = 
| /i, 777,^)1/2, to/ 2 ). Here, fj = Sj + ij is the total spin 
of atom j and mj its projection onto the -B-field axis. 

In the following we consider an elastic collision, where 
only the entrance channel with spin configuration \ao) is 
open and all other coupled channels are closed. 

Within the TC approximation of the scattering pro- 
cess one projects the full MC Hilbert space onto two 
subspaces, the one of the closed channels (with projec- 
tion operator Q) and the one of the open entrance chan- 
nel (with projection operator P) The resulting TC 
Schrodingcr equation reads 

(U P - E)\y P )+W\y Q ) = (4) 

$Q-E)\9 Q )+Wt\9 P ) = 0, (5) 

with Hp = PHP, H Q = QHQ, W = PHQ, |* P ) = P|*) 
and =Q|*). 

An MFR occurs, if the energy E of the system is close 
to the eigenenergy E of a bound state |$&) of the closed- 
channel subspace. Following the solution in [l3[ we as- 
sume that the closed-channel wave function is simply a 
multiple A of the bound state |<&j,), i.e. \9q) = A|<&&). 
This is equivalent to the usual one-pole approximation. 
Equation (j4]) may be solved via the Greens operator 
G P = (E + ie- Hp)^ 1 with e -> 0+ 

The general solution thus reads 

\*p) = C\<f> ice ) + AG P W\<f> b ) , (6) 
l*Q> - A\^ b ) (7) 
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where C is a normalization constant and |3> re g) is the 
regular solution of Hp I*) = E\^). 

Following the procedure given in [l3T ] one arrives at the 
closed-channel admixture 



.4 



-C 



sin 6 r , 



(8) 



with the normalization constant C = C'/ cos5 rcs . The 
Greens operator Gp explicitly given in [H yields a long 
range behavior of the open channel 



Selection rules allow only transitions from the s-wave 
scattering function to a final state with J = 1. Due to 
the orthogonality of the molecular basis, only one molec- 
ular channel i\) Xi (R) with the same spin state as the final 
state has to be considered. 

The solutions (|6l7p of the TC approximation yield to- 
gether with the behavior of the closed channel admixture 
in ([5]) a squared dipole transition moment 



I TC (B) = |(* f |D|*)| 



| C\ 2 | Ci cos 5 ICS - C 2 sin 5 T , 



(14) 



Vp(R)\ 



R — >oo 



= c 



2fi 
irh 2 k 



sm(kR + S hg + 5 ICS ) . (9) wh ere C 2 = v^/ttT ((* f |DG P W|$ 6 ) + (* f |D|$ 6 » and 



The total phase shift 8 — <5bg + S Ies results from the 
background phase shift (5b g of the regular solution |$ re g), 
which is connected to the background scattering length 
dbg = — linifc^o tan^bg/fc, and a contribution 5 res due to 
the resonant coupling to the bound state. The resonant 
phase shift has the functional form 



tan<5 rns = — - 



r/2 

E — Er, 



(10) 



where E R = E + ($ b |W t GpW|<l>f ) ) lies close to the en- 
ergy of the bound state E (B). The width of the reso- 
nance is given by T = 2tt\ (${,|W|$ rcg ) | 2 . 

One assumes that Er depends approximately linearly 
on the magnetic field, i.e. for E — > there is some a 
and Br such that Er = a(B — Br). This yields with 
AB = T{2ka\ >s (i)^ 1 the well known relation [lj| 



a bg 1 



AB 

B-Br 



which allows to determine the resonant phase shift 



tan <5 r , 



fccibgAi? 

Br — B 



(11) 



(12) 



from experimentally accessible quantities. 

Equipped with the MC solution a convenient way to 
calculate transition rates to molecular bound states is to 
transform the scattering wave function of Eq. §5§ to the 
molecular basis \\) — \S, Ms) \rrii 1 , rrn 2 ) where S and Ms 
are the quantum numbers of the total electronic spin and 
its projection along the magnetic field and , m l2 are 
the nuclear spin projections of the individual atoms. 

Within the dipole approximation with electronic 
dipole transition moment D(R) the free-bound tran- 
sition rate T(B) to the final molecular state |^f) = 
^ r "j^ Yj I (@, $)|xf) with vibrational quantum number v 
and rotational quantum number J is then proportional 
to the squared dipole transition moment [1 51 ] 



Imc(B) = 



oc 

I 



V u (R)D(R)tp Xl (R)dR 



(13) 



C\ = ("ft |D|$ rcg ) do not vary with the magnetic field 
B. The prefactor C may vary with a sc . However, in 
the following we consider energy-normalized scattering 
solutions with C = 1. Introducing (3 = \C\ 2 (C 2 + Cf) 
and tan<5o = C\jCi one can further simplify Eq. (I14|) to 



Itc(B) =/3sin 2 (5« 



So) 



(15) 



Due to the low-energy behavior of the regular solution 
r cx k and C\ cx Vk holds for k — > ■ Therefore one 
can associate a finite length 



lim 

k^O 



tan 8q 



(16) 



with the phase shift <5o. 

Let us point out some differences to the previously de- 
rived result for the dipole transition moment in [5j. In 
the notation of the current work Eq. (8) in Q gives 

Jtc(B) =X|l + CitanJ rcs + C 2 sin ( 5 rcs | 2 (17) 



with K = |(#f|D|$ r 
and C*2 = 



d = (* f |D|$ irr )/(* f |D|$ 



rcg/, 



x /27^r(* f |D|$ b )/(* f |D|$ reg ). The most 



obvious difference to Eq. (|14p is the dependence on three 
parameters K, C\ and Ci and not just two. This is a re- 
sult of an inconsistent normalization of open and closed 
channels in [5[. The open channel was not energy nor- 
malized and leads thus to a term proportional to tan i5 ros . 
Furthermore, the open channel was described as a pure 
sum of regular and irregular solution. This may, however, 
only be done for interatomic distances, were the coupling 
to the closed channels induced by the exchange energy is 
negligible. A fit of Eq. (jXTJ) to a full MC calculation 
seemed to be nevertheless possible which might however 
be the result of the freedom of three fitting parameters. 

The universal dependence of the PA transition rates 
on just two parameters in Eq. (|15[) reveals an important 
physical aspect. The enhancement of the transition rate 
is directly connected to the position of vanishing transi- 
tion rate. On the one hand the transition rate is van- 
ishing at a scattering length ay™ = a\ 
corresponding magnetic field 



a e or at a 



B„ 



Br 



AB 



(18) 
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On the other hand the point of vanishing transition rate 
is connected to the enhancement ratio 



1 



1 



bg 



1 S k 2 a 2 e 



-B f 



ka hR AB 



(19) 



of the maximum transition rate r max and the background 
rate Tb g in the presence of an off-resonant magnetic field 
(i.e. where <5 rcs = 0). 

In order to verify the TC description of the PA process, 
we consider the exemplary case of an elastic collision of 
6 Li- 87 Rb ( 6 Li is atom 1, 87 Rb is atom 2) in the initial 
atomic basis state |ao) = 1 1/2 , 1/2) 1 1 , 1). For an energy 
50 Hz above the threshold of the entrance channel which 
is well in the s-wave scattering regime the MC solution 
was calculated for different magnetic fields B in [l6[ . For 
B < 1500 G two s-wave resonances occur, a broad one at 
B = 1066, 917 G which was also observed experimentally 
and a narrow one at B = 1282.576 G. The depen- 
dence of the scattering length a sc on the magnetic field 
strength is shown in Fig. [1] 
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FIG. 1: (Color online) Scattering length a BC as a function of the 
external magnetic field value B for 6 Li- 87 Rb scattering at E = 
50 Hz (dots). A fit according to Eq. H20H is depicted by the solid 
(red) line. The value of a]~, g (B) is shown by the dashed (orange) 
line. All fitting parameters are shown in the plot. 

Assuming that the two resonances are sufficiently sepa- 
rated in order to describe the process by two independent 
resonances one may generalize Eq. pip to 



a sc (B) = Ob g I 1 



AB 1 



AB 2 



B — Bri B — Br2 



(20) 



Additionally, we account for effects beyond the one-pole 
approximation by allowing ab g to vary slowly with B as 
abg(£>) = ao + ai • B + a 2 ■ B 2 . With this quadratic 
expansion, a fit according to Eq. ([2D)) excellently reflects 
the MC behavior as shown in Fig. [TJ 

We consider the exemplary case of a dipole transitions 
of the scattering state to the absolute vibrational ground 
state of the electronic singlet configuration AT 1 E + and the 
triplet configuration a 3 E + . These transitions take place 
at internuclear distances where the coupling between all 
atomic channels is strong such that any deficiency of the 
TC description would become obvious. The MC rate 
was calculated in [l6j for an electronic dipole moment 
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FIG. 2: (Color online) Squared dipole transition moment as 
a function of the external magnetic field for 6 Li- 87 Rb scatter- 
ing at E = 50 Hz (dots). Transitions to the molecular singlet 
state |0, 0) |0, 3/2) (a) and the triplet states |1, — 1> [1, 3/2) (b) and 
1, 1)|1, —1/2) (c) are considered. In each graph a fit according to 
Eq.s 12 1 [ i and i !22[ i is performed and the resulting fitting parameters 
for Eq. (12 1 1 are shown. The insets focus respectively on a region 
B R2 - 2AB 2 < B < B R2 + 2AB 2 . 



in the linear approximation D(R) — Dq + D± ■ R, where 
D could be neglected. In the following we use Do = 
and Di = Eh/ao. The magnetic field dependence of the 
according dipole transition moment to ground states in 
different spin configurations is shown in Fig. O 

We fit the MC behavior by assuming again that the two 
resonances are sufficiently separated, such that one can 
add the transition amplitudes Aip c (B) — Co cosmos ) — 
Cj-sin^) of both resonances j — 1,2 independently 
and take the absolute square of the sum 



4 1 c 2) (S) = |^(i?) + /4^( J B)| 



id) 



( 2 ), 



(21) 



in order to determine the dipole transition moment. The 

resonant phase shifts fancies = kab g {B)ABj/(Bnj — B) 
for j = 1,2 are in analogy to Eq. (fT2"|) associated to the 
resonant coupling to two different closed-channel bound 
states. Since Co does not depend on the resonant molec- 
ular bound state it is independent of the magnetic 
field B. Hence, for describing a transition process to a 
specific molecular state for two well separated resonances, 
one needs only three independent parameters. 
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TABLE I: Enhancement of the dipole transition rate be- 
tween off-resonant and resonant magnetic field in MC and 
TC description for transitions to all eight possible spin states 
\S, Ms)\m,i 1 ,rrii 2 ) at the two resonances at Br\ = 1066, 917 G 
and B R2 = 1282.576 G. 



Molecular 


Resonance 1 


Resonance 2 


state 


7-max ; 7-oft (U „W\ 
J MC /- l MC ( K a e I 


-2 rmax ; joS f I, „ ( 2 ) \ -2 



|0,0)|0,|) 


5.81 


10 6 


5.83 


10 6 


6.85 


10 12 


1.39 


10 13 


|0,0>|1,|) 


2.62 


10 8 


3.49 


10 8 


6.57 


10 6 


1.27 


10 7 


|x,-i)|i,|) 


1.80 


10 5 


1.75 


10 5 


5.73 


10 7 


5.41 


10 7 


[1, 0>|0, §> 


8.86 


10 6 


8.69 


10 6 


1.63 


10 13 


5.93 


10 13 


|1,0>|1,|) 


5.04 


10 9 


2.14 


10 9 


2.15 


10 8 


4.71 


10 7 


|1, 1)1-1,1) 5.77 


10 6 


8.41 


10 6 


7.76 


10 12 


5.21 


10 13 


|1,1>|0,|) 


1.88 


10 7 


2.13 


10 7 


1.32 


10 13 


6.57 


10 13 


|X, 1>|1, — i> 


1.09 


10 10 


1.18 


10 12 


3.66 


10 7 


4.41 


10 7 



Additionally a fit to the behavior 



j) 

res 



(22) 



for j = 1,2 is performed which neglects respectively one 
resonance. Again one can connect the phase shifts 5^ 

to the corresponding lengths ai via Eq. (|16p which al- 
lows also to validate the applicability of Eq. (fl~9]) for both 
resonances separately. 

Considering Fig. [2] one finds that for all shown transi- 
tions the TC approximation for two well separated res- 
onances excellently describes the magnetic-field depen- 
dence of the MC transition rate. The behavior of Eq. (j22|) 
reproduces each of the MC resonances j = 1, 2 in an inter- 
val of several ABj around the resonant magnetic fields 
Brj. With only one parameter more than for the de- 
scription of a single resonance, both resonances can be 
well described by Eq. (f2Tj) . This is remarkable, since the 
behavior of the transition rates to all four different spin 
configurations are quite different. Nevertheless all fea- 
tures are reproduced by just three parameters. 



The good description of the MC transition rates by 
Eq. (|22|) suggests that Eq. (fl~9"|) indeed reflects the de- 
pendence of the PA enhancement on the position of van- 
ishing transition rate. In Eq. (|19p the enhancement was 
defined relative to the transition rate at the background 
scattering length. This point is reached, however, only 
at infinite detuning from the resonant magnetic field Br. 
In order to nevertheless verify the validity of Eq. (|19p we 



relate the maximal transition rate /{JcT a ^ eacn resonance 



separately to the transition rate Iyf c = Imc(800G) far 
away from both resonances. We do not choose a mag- 
netic field with larger detuning to avoid effects of other 
resonant molecular bound states. In Tab. [J the ratio 
Imc/Imc i s compared for both resonances to the predic- 
tion of Eq. (I19|) for transitions to the vibrational ground 
states of all eight possible spin configurations in the 
molecular basis. One finds that the order of magnitude 
generally agrees excellently. Only for few transitions such 
as the one to the molecular states |1, 1)|1, — 1/2) at the 
first (broad) resonance the orders of magnitude differ sig- 
nificantly. A view on Fig. [2] (c) reveals that this is not 
related to a break down of Eq. (fl~9"l) . but that the ab- 
sence of a vanishing transition rate leads to a comparably 
slow degradation of the transition rate such that I^c ^ s 
not a good representation for the background transition 
rate. On the other hand, for transitions for which the 
background transition rate is quickly approached when 
detuning from the resonance, the two estimates of the 
enhancement agree even to the first significant digit (see 
the third row in Tab. |l] and the corresponding Fig.[2](b)). 
This and the results above demonstrate that the TC ap- 
proximation provides an excellent basis to understand PA 
processes in the presence of an external magnetic field in- 
ducing an MFR. 
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